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ABSTRACT 

Hubble Space Telescope (HST) WFPC2 archival images obtained years prior to the explosion of the 
Type Il-plateau supernova (SN) 2001du in NGC 1365 previously have been analyzed to isolate the 
progenitor star. The SN site was located using precise astrometry applied to the HST images, with 
significant uncertainty, leaving three possible progenitor candidates. Images of the fading SN have 
recently become publicly available in the HST archive, allowing us to pinpoint the SN's exact location 
on the pre-explosion images. We show that the SN occurred in very close proximity to one of the blue 
candidate stars, but we argue that this star is not the actual progenitor. Instead, the progenitor was not 
detected on the pre-SN images, and we constrain the progenitor's mass to be less than 13+4 M Q . This 
is consistent with previous constraints on the progenitor masses of other Type Il-plateau supernovae 
(SNe), suggesting that such SNe arise from the iron core collapse of massive stars at the lower extreme 
of the possible mass range. 

Subject headings: supernovae: general — supernovae: individual (SN 2001du) — stars: massive — stars: 
evolution — stars: variables: other — galaxies: individual (NGC 1365) 



1. INTRODUCTION 

Supernova 2001du was visually discovered by R. Evans 
(2001) on August 24.7 (UT dates are used throughout 
this paper) at the western end of the bar in the spiral 
galaxy NGC 1365 (about 90" W and 10" S of the nucleus). 
Smartt, Kilkenny, & Meikle (2001) classified the SN as 
Type II; Wang et al. (2001), from a spectrum obtained 
on August 30 UT, confirmed the SN as Type Il-plateau 
(II-P) before optical maximum. NGC 1365 has also been 
host to SN 1957C (of unknown type) and the Type Ic SN 
1983V. Using Hubble Space Telescope (HST) observations 
of Cepheids, Silbermann et al. (1999) determined the dis- 
tance modulus of NGC 1365 to be /i = 31.3±0.2 mag, 
corresponding to a distance of about 18.2 Mpc. 

Van Dyk, Li, & Filippenko (2003) recently measured 
the absolute position of SN 2001du and, applying an inde- 
pendent astrometric grid to pre-SN HST archival WFPC2 
images, isolated three possible candidate progenitor stars 
within the uncertainty of the measured SN position: one 
reddish star with a detection in both the F555W and 
F814W bands, and two stars with only F555W detections, 
implying bluer colors for these latter two stars. Also, 
multi-band WFPC2 images of the SN were obtained for 
HST program GO-9041 (PI: S. Smartt) in 2001 Novem- 
ber. These images are now in the public domain, and 
in this paper we analyze them and compare them to the 
pre-SN images to further constrain the nature of the SN 
progenitor. (Smartt et al. 2003 have also recently analyzed 
their images and, in a somewhat similar manner, come to 
the same basic conclusion that we reach here.) 



2. ANALYSIS 

Pre-SN images, each of total exposure time 300 s in 
the bands F336W, F555W, and F814W, were obtained by 
GTO-5222 on 1995 January 15. Late-time multi-band im- 
ages of the SN were obtained by GO-9041 on 2001 Novem- 
ber 26. These include deep images in F336W, F439W, 
F555W, F675W, and F814W, as well as shorter exposures 
in F336W, F555W, and F814W. The available HST data 
are summarized in Table 1, with total exposure times given 
for each band. 

In Van Dyk et al. (2003) we adopted the position of 
the SN, measured from an image obtained with the YALO 
1-m telescope at CTIO (and available on the Internet), 
as a = 3 ft 33 m 29!15, 6 = -36°08'31"5 (J2000), with to- 
tal uncertainty ±0"9. The SN site is located on the WF3 
chip. In Figure 1 we reproduce Figure 21 from Van Dyk 
et al., showing the SN environment in the F555W band. 
In Figure 2 we show the late-time F555W image (the SN 
is nearly centered on the PC chip for all the GO-9041 ex- 
posures) . As one can see, the actual SN position is within 
the 0'.'9 error circle, which provides additional confidence 
in the astrometric method used by Van Dyk et al. when 
isolating the various SN sites 2 . 

We initially used the routine HSTphot (Dolphin 
2000a, b) to obtain photometry of the objects imaged 
on the PC chip. HSTphot automatically accounts 
for WFPC2 point-spread function (PSF) variations and 
charge-transfer effects across the chips, zeropoints, aper- 
ture corrections, etc., and can return as output magni- 



1 BASED ON OBSERVATIONS MADE WITH THE NASA/ESA HUBBLE SPACE TELESCOPE, OBTAINED FROM THE DATA 
ARCHIVE OF THE SPACE TELESCOPE SCIENCE INSTITUTE, WHICH IS OPERATED BY THE ASSOCIATION OF UNIVERSITIES 
FOR RESEARCH IN ASTRONOMY, INC., UNDER NASA CONTRACT NAS 5-26555. 

2 Applying the same astrometric method used by Van Dyk et al. (2003) to the mosaic of the deep F555W exposure, we measure a position 
for SN 2001du to be a = 3 fe 33 m 29?14, <5 = -36°08'31('3 (J2000), ±0'.'2 in each coordinate. 
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tudes transformed from the HST flight system to stan- 
dard Johnson-Cousins bands. Van Dyk et al. (2003) dis- 
cuss the similarity in the results from HSTphot, compared 
with those from the packages DoPHOT (Schechter, Mateo, 
& Saha 1993) and DAOPHOT / ALLSTAR (Stetson 1987, 
1992). 3 

The SN profile is saturated in the longer F555W (V), 
F675W (R), and F814W (J) exposures. Thus, no bright- 
ness measurement is possible. The SN is not saturated 
in the shorter F555W and F814W exposures, the long 
F336W (U) and F439W (B) exposures, and also the short 
F336W exposure. However, the results from HSTphot for 
the F336W exposures varied, depending on how these ex- 
posures were combined with the other bands as input into 
the routine, and the results from the long and short expo- 
sures for this band were also discrepant with each other 
(the results ranged from about mF336W = 18-9 to 18.7 mag 
for different runs of HSTphot). We do not know why this 
occurred, but we speculate that it possibly may be due to 
the low signal-to-noise ratio for both exposure sets. For 
F336W we then applied to both the long and short expo- 
sures DAOPHOT/ ALLSTAR within IRAF 4 as a check, us- 
ing a TinyTim PSF (Krist 1995), the zeropoint from Holtz- 
man et al. (1995), and aperture corrections (^0.00 mag) 
derived from artificial star tests. (Checks for the F439W 
and short F555W and F814W exposures with DAOPHOT 
show good agreement with the HSTphot results.) 

We obtained to F3 36W = 18.96±0.10 and 18.87±0.04 mag 
for the short and long exposures, respectively. We then av- 
eraged these DAOPHOT magnitudes with' the HSTphot 
magnitudes at F336W, and obtained mF336W = 18.86 
mag. The SN was quite red by late November 2001, with 
WF336W — "7-F439W = 1.82 mag. Whether using the trans- 
formations from flight system to Johnson-Cousins magni- 
tudes from Holtzman et al. (1995), Dolphin (2000b), or 
via SYNPHOT applied to the Bruzual Spectral Synthetic 
Atlas (see Van Dyk, Filippcnko, & Li 2002), this color is 
too red to obtain a reliable transformation from TOF336W 
to U. However, we do apply a transformation, with uncer- 
tainty ^0.2 mag. The transformations for all other bands 
are from Dolphin (2000b), applied internally within HST- 
phot. The final results for the photometry of SN 2001du 
in the HST images are presented in Table 2. 

In Figure 3 we show preliminary, unpublished BVI light 
curves for SN 200 ldu, based on data from the YALO 1- 
m telescope (N. Suntzeff & K. Krisciunas, 2003, private 
communication), and we include the UBVI photometry 
from the HST data. For comparison we show the B and 
V light curves for the plateau-type SN (II-P) 1999em in 
NGC 1637 (Leonard et al. 2002a). It is clearly evident that 
SN 2001du is indeed a SN II-P as well. If this SN were a 
linear- type event (ILL), it should have markedly declined 
in brightness by the time of the HST observations (~94 
days since discovery). 

3. DISCUSSION 



In Van Dyk et al. (2003) we had isolated three faint 
objects, A-C, detected and determined to be stellar by 
HSTphot, within the error circle around the SN position. 
Stars A and B, to the south and west, respectively, are 
both blue, with ttt-fsssw (~V) = 24.30±0.21 mag and 
25.02±0.32 mag, respectively, and no detections at F814W 
(~/). (Blue, here, is relative, and includes stars with 
roughly early-G spectral type or earlier, based on the de- 
tection limits; see below.) Star C, to the east, is relatively 
red, with V = 24.44±0.23 and V - I = 1.03±0.30 mag. 
We had assumed that, since SN 2001du is of Type II-P, 
the most plausible of the three candidates is the redder 
star. However, that was argued to be still unlikely, based 
on the assumed luminosity and color of Star C. 

We have used the late-time SN images, specifically the V 
image, and attempted to match them up with the pre-SN 
images. This was done by using the task geotran in IRAF 
to transform both the late-time SN images to the pre-SN 
images and vice versa. We can therefore locate the position 
of the SN on the pre-SN images as being 0.70±0.15 WF 
pixel (0'.'07) northeast of the position of Star B. It would 
be tempting to assign Star B, then, as the progenitor of 
SN 2001du. However, (1) the offset of the SN position 
from Star B (which corresponds spatially to 6.2 pc at the 
distance of NGC 1365) is significantly large, and the accu- 
racy of the measurement is appreciably high; (2) SNc II-P 
are thought to arise from the explosions of red supergiants; 
and (3) no evidence exists that SN 2001du is peculiar in 
any way and that it arose from a bluer, more compact 
star like the progenitor of SN 1987A (E. Cappellaro, 2002, 
private communication). Thus, it is more likely, as we 
previously had suspected (Van Dyk et al. 2003), that the 
progenitor simply is not detected in the pre-SN images. 

We can use the detection limits of both the F555W and 
F814W pre-SN images to constrain the nature of the pro- 
genitor. We do this by inferring limits on the luminosity 
and, therefore, the mass from the limits on the star's ap- 
parent magnitude. These HSTphot detection limits in the 
SN environment are mF555\v > 25.19 and mF8i4W > 24.25 
mag. Before converting into absolute magnitude, we must 
first assess the amount of extinction to the SN. In Van Dyk 
et al. (2003) we attempted to use stars in the environment 
to estimate the reddening. The Galactic extinction to- 
ward NGC 1365 is A v = 0.07 mag (Schlegel, Finkbeincr, 
& Davis 1998; and NED 5 ). We had found that stars in 
the larger SN environment implied larger amounts of ex- 
tinction, Ay rts 2.5 mag. (Unfortunately, although the 
post-SN images are deeper than the pre-SN ones, the SN 
is far too bright, and its image is spread over too many 
surrounding pixels, to perform a more thorough study of 
the environment.) 

Using the light curve in Figure 3 we can actually derive 
more specific information about the possible reddening to- 
ward the SN from the SN itself (which is information that 
was not available to us previously). The HST photometry 
of the SN on November 26 is about 94 days after discov- 
ery, i.e., near the end of the light-curve plateau, where we 



3 In fact, here we update the private communication by D. C. Leonard 2002, regarding the very good agreement between the results of 
HSTphot photometry on images of NGC 3351, compared with those obtained by Graham et al. (1997) using DAOPHOT; Leonard now finds 
5V = 0.015±0.022 and 51 = 0.036±0.018 mag. 

4 IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 

5 NED is the NASA/IPAC Extragalactic Database, http://ncdwww.ipac.caltech.edu. 
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can still estimate the SN's color temperature with reason- 
able accuracy (e.g., Hamuy et al. 2001). These colors are 
B — V = 1.42 and V — I = 0.91 mag. For comparison, the 
colors of the very well-studied Type II-P SN 1999em at the 
same age (about 95 days after discovery) are B — V = 1.47 
and V — I = 0.89 mag (Leonard et al. 2002a). Leonard 
et al. and Hamuy et al. (2001) conclude that SN 1999em 
is reddened by E(B — V) = 0.10 mag. Since B - V is 
affected by line blanketing in the B band (e.g., Hamuy et 
al. 2001), the small difference in B- V (SN 1999cm is 0.05 
mag redder than SN 2001du) could be due to differences in 
metallicity between the two SNe, which might affect this 
color. The marginally redder V — I color of SN 2001du 
(by 0.02 mag) may imply that this SN is slightly more 
reddened than SN 1999em. However, we consider these 
differences in color between the two SNe to be small and 
likely within the uncertainties in the colors, so we adopt 
E(B - V) « 0.1 mag for SN 2001du, as found for SN 
1999em. 

If we assume the Cardelli, Clayton, & Mathis (1989) 
reddening law, then Ay w 0.3 mag and Ai « 0.2 mag 
for SN 2001du. For the distance modulus \x — 31.3 mag, 
My > —6.4 and Mj > —7.0 mag for the progenitor star. 
In a similar vein to the estimation of the progenitor mass 
limit for SN II-P 1999em by Smartt et al. (2002) and SN 
II-P 1999gi by Smartt et al. (2001), the absolute magni- 
tude limits can be converted to the luminosity of the likely 
supergiant progenitor, assuming the full range of possible 
stellar surface temperatures. This is accomplished using 
a similar relation to equation (1) in Smartt et al. (2001, 
2002), where we also have assumed the temperatures and 
stellar bolometric corrections for supergiants from Drilling 
& Landolt (2000) and have computed, using SYNPHOT 
within IRAF and the Bruzual Spectral Atlas, the conver- 
sions of tof555w and mF8i4W to Johnson V for the range 
in spectral type. The luminosity limits for both F555W 
(V) and F814W (J) are shown in Figure 4. Stars with lu- 
minosities brighter than these limits (shown by the heavy 
lines) should have been detected in the pre-SN images. 

For comparison we have plotted (also in a similar fashion 
to Smartt et al.) the model stellar evolutionary tracks for 
a range of masses from Lejeune & Schaerer (2001), assum- 
ing enhanced mass loss for the most massive stars and that 
a metallicity Z = 0.04 is appropriate for the SN environ- 
ment. Indication for a metallicity somewhat greater than 
solar in this sector of NGC 1365 comes from the abundance 
study by Roy & Walsh (1997) of H II regions in the galaxy. 
The closest H II regions to the SN site are their Nos. 20 
and 21, where at this radial distance from the nucleus, 
the O/H abundance is ~9 dex or, in fact, slightly greater, 
based on model assumptions (the solar O /H abundance is 
8.8 dex; Grevesse & Sauval 1998). 

The luminosity limits at each band in Figure 4 show that 
the pre-SN images are not sensitive to the bluest high-mass 
supergiants, so we cannot be absolutely sure that such a 
star is not the progenitor. If we exclude this possibility, 
the V band sets an upper limit to the progenitor mass 
of ~ 13 Mq, although this ranges upward to as high as 
~ 20 Mq, depending on the effective temperature. The 
I band is more sensitive to less-massive red supergiants, 
ranging downward to <~ 9 Mq (possibly even less) at the 
lowest effective temperatures. These results from Figure 4 



are roughly consistent with what we showed in Figure 22 
of Van Dyk et al. (2003). We therefore estimate that the 
upper mass limit for the SN 2001du progenitor is 13+4 

This limit is consistent with the limit Leonard et 
al. (2002b) find for SN II-P 1999gi (15±| Mq) and the 
- 17 Mq limit one finds for SN II-P 1999cm, if the dis- 
tance to the SN 1999em host is assumed to be the Cepheid 
distance, 11. 5±0. 9 Mpc (Leonard et al. 2003). Low masses 
for the progenitors of SNe II-P, relative to that indicated 
for SNe II-L (linear) and SNe Iln (narrow), are also in- 
ferred from their optical properties (e.g., Schlegel 1996; 
Filippcnko 1997; Chugai & Utrobin 2000) and suggested 
by the inferred paucity of circumstellar matter from radio 
and X-ray observations (e.g., Weiler et al. 1989; Schlegel 
2001; Pooley et al. 2002). 

4. CONCLUSIONS 

Comparing HST archival WFPC2 images obtained of 
the Type II-P SN 2001du in NGC 1365 at late times with 
images obtained before explosion, we have located the po- 
sition of the SN progenitor star to within 1 WF chip pixel 
of the candidate Star B identified by Van Dyk et al. (2003). 
This star, however, is more than likely not the progeni- 
tor. We use the detection limits of the pre-SN F555W and 
F814W images to place constraints on the progenitor's na- 
ture; specifically, M <; 13±4 Mq. (Smartt et al. 2003 find 
a similar mass limit of M = 15 Mq.) This limit is con- 
sistent with what has been found and inferred for other 
SNe II-P. 

The proximity of the SN to the blue Star B implies that 
the progenitor may have been part of a small stellar associ- 
ation (the separation between the progenitor and Star B is 
likely too large for the blue star to be a binary companion). 
The faint emission underlying Stars A, B, and C identi- 
fied by Van Dyk et al. (2003) suggests that an association 
complex may lie in the SN environment. It would be inter- 
esting to image the SN again in multiple bands with HST 
at high spatial resolution, presumably with the Advanced 
Camera for Surveys (ACS), when the SN has significantly 
dimmed. The stars in the environment of the fading SN 
would then become more obvious, and the age and mass 
of the progenitor star might be further constrained, based 
on the characteristics of its surviving neighbors. 

This research made use of the NASA/IPAC Extra- 
galactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy, under contract with NASA, and the LEDA database 
(http://leda.univ-lyon.fr). The work of A.V.F.'s group 
at UC Berkeley is supported by NASA grants AR-8754, 
AR-9529, and GO-8602 from the Space Telescope Sci- 
ence Institute, which is operated by AURA, Inc., under 
NASA contract NAS5-26555. We thank D. C. Leonard and 
M. Hamuy for useful discussions. We also thank the ref- 
eree, K. Krisciunas, for helpful comments and for provid- 
ing unpublished ground-based photometry of SN 2001du. 
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Table 1 

Summary of HST Observations 



Date 


Filter 


Exp. Time 


HST 


SN 


(UT) 




00 


Program 


Saturated? 


1995 Jan 15 


F336W 


300 a 


GTO-5222 






F555W 


300 








F814W 


300 






2001 Nov 26 


F336W 


100 


GO-9041 


N 






460 




N 




F439W 


400 




N 




F555W 


40 




N 






600 




Y 




F675W 


400 




Y 




F814W 


40 




N 






600 




Y 



a The exposure times for all three bands in this 1995 Jan 15 dataset were erroneously listed in Van Dyk et al. (2003) 
as 100 s. The actual total exposure time of 300 s results from the combination of a trio of 100-s exposures. Note that 
HSTphot automatically reads the correct total exposure time from the images during processing, so the reported results 
are unaffected by this error. 
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Table 2 

Photometry of the SN 2001du ffST Images 



u 


B 


V 


R 


I 


[TOF336W] 


[WF439W] 


[™F555w] 


[WF675W] 


[TOF814W] 


18.6(0.2) a 


16.72(0.06) 


15.30(0.00) 




14.39(0.00) 


[18.86] b 


[17.04] 


[15.32] 




[14.43] 



a The value in parentheses for each entry in this row is the uncertainty in the transformed Johnson-Cousins 
magnitude. 

b The value in brackets for each entry in this row is the adopted flight-system magnitude, derived from 
HSTphot, DAOPHOT, or both. 



Fig. 1. — The site of SN 2001du in the archival 300-s F555W image from 1995 January 15. The error circle has radius 
0'.'9. The three stars, A-C, within the circle are indicated with tickmarks. This is a reproduction of Fig. 21 from Van 
Dyk et al. (2003). 



Fig. 2. — SN 2001du, as seen in the 600-s F555W image from 2001 November 26, shown to the same scale and orientation 
as in Figure 1. The error circle from Figure 1 is also shown. The supernova is saturated in this image. 
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Fig. 3. — The light curves of SN 2001du. The open pentagons, open squares, and open triangles are BVI measurements 
from images obtained using the YALO 1-m telescope (N. Suntzeff & K. Krisciunas, 2003, private communication). The 
filled symbols are the measurements from the HST images. Also shown for comparison are the measurements in B (five- 
pointed stars) and V (four-pointed stars) for the Type II- P SN 1999em from Leonard et al. (2002a), adjusted by ~1 mag 
to match the V peak for SN 2001du. 
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Fig. 4. — Hcrtzsprung-Russcll diagram showing model stellar evolutionary tracks (alternating long-dashed lines and short- 
dashed lines) for a range of masses from Lejeune & Schacrer (2001), with enhanced mass loss for the most massive stars 
and a metallicity Z = 0.04. Also shown are the luminosity limits (heavy solid lines) for both the F555W (V) and F814W 
(I) bands. Stars with luminosities brighter than these limits should have been detected in the pre-SN HST images. 
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